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Periodic class II methanol masers in G9.62+0.20E 
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ABSTRACT 

We present the light curves of the 6.7 and 12.2 GHz methanol masers in the star 
forming region G9.62+0.20E for a time span of more than 2600 days. The earlier 
reported period of 244 days is confirmed. The results of monitoring the 107 GHz 
methanol maser for two flares are also presented. The results show that flaring occurs 
in all three masing transitions. It is shown that the average flare profiles of the three 
masing transitions are similar. The 12.2 GHz masers are the most variable of the three 
masers with the largest relative amplitude having a value of 2.4. The flux densities 
for the different masing transitions are found to return to the same level during the 
low phase of the masers, suggesting that the source of the periodic flaring is situated 
outside the masing region, and that the physical conditions in the masing region 
are relatively stable. On the basis of the shape of the light curve we excluded stellar 
pulsations as the underlying mechanism for the periodicity. It is argued that a colliding 
wind binary can account for the observed periodicity and provide a mechanism to 
qualitatively explain periodicity in the seed photon flux and/or the pumping radiation 
field. It is also argued that the dust cooling time is too short to explain the decay time 
of about 100 days of the maser flare. A further analysis has shown that for the intervals 
from days 48 to 66 and from days 67 to 135 the decay of the maser light curve can be 
interpreted as due to the recombination of a thermal hydrogen plasma with densities 
of approximately 1.6 x 10^ cm~'^ and 6.0 x lO'^ cm^"^ respectively. 
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1 INTRODUCTION 

Although class II methanol masers are now generally ac- 
cepted to be exclu sivel y associated with massive s tar form- 
ing regions (see eg. lEilingsen . 20061 ; IXu"et aLlbOOSl ). it is not 
yet clear how much can be learned about the star formation 
environment from the masers. Underlying this uncertainty 
is the fact that it is difficult to determine where in the cir- 
cumstellar environment the masers operate based solely on 
the spatial distributi on and velocity structure of the masers 
l|Beuther et al.ll2002l ). In the past very little attention has 
been given to the time domain aspect of the masers. Var- 
ious a uthors (see eg. MacLeod ct a f] 1 19931; ICaswell et al.l 



19951; iMoscadeUi fc Catarzii |1996|; iCocdhart ct al.l |2002|; 



Niezurawska et al.ll2002 ; ,Goedhart et al.,.200a . .20041 . l2005al ) 
have studied the variability of the 6.7 and 12.2 GHz masers 
over various periods of time and have found that variabil- 
ity is a common feature amongst these masers. The most 
systematic study on the variability of 6 .7 GHz masers to 
date was done by iGoedhart et al. I l|2004l ) using the Harte- 
beesthoek Radio Astronomy Observatory (HartRAO) 26-m 
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telescope, revealing a wide variety of time-dependent be- 
haviour in the masers. Six out of their list of 54 sources were 
identi fied as periodic, with periods ranging from 133 to 504 
days l|Goedhart et al.ll2007l ). Of particular interest here are 
the periodic masers in the star forming region G9. 62-1-0. 20E 
which sho w repeated flaring activity with a period of about 
244 days (|Goedhart et al.ll2003l ; iGavlard fc Goedhart|[2007h 

The star forming region G9. 62-1-0. 20 has been t he object 
of study by a number of authors. iGarav et al.l (|l993l ) mapped 
this region with the VLA at 1.5, 4.9 and 15.0 GHz and iden- 
tified five extended, compact and ultracompact (UC) Hll 
regions, labelled A - E. A striking feature in this star form- 
ing complex is the near perfect alignment of various tracers 
(masers, hot molecular clumps and UC H II regions) of star 
forming activity extending between compone nts C and D, 
with c omponent E l y ing b etween C and D (see lHofner et all 
ll994h . [Hofiier et al.l l|l994h interpreted the overall morphol- 
ogy of this star forming region as suggestive of induced 
star formation progressing from the most evolved H II region 
(component A) with the most recently formed stars located 
in the linear structure. The methanol masers in the star 
forming complex are associated with components D and E 
and are in the literature also referred to as G9. 619-1-0. 193 
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and G9.621+0.196 respectively l|Pliillips et all 1 19981 ). For 
component E, the masers have a maximum linear extent 
of about 566 AU (2.75 milliparscc) based on the 12.2 GHz 
observations of .Goedhart et a l. (2005b) and using a distance 
of 5.15 kpc (A. Sanna et al., in preparation). 

Component E shows all the signs of a very early phase 
of massive star formation. According to iHofner et al.1 ()l996l l 
the continuum emission follows a power law with a spec- 
tral index of 1.1 ± 0.3 between 2 cm and 2.7 mm, which 
can be interpreted to be due to either an ionized spheri- 
cal stellar wind or a spherical homogeneous UC Hll region 
with an excess of dust emission at 2.7 mm. If interpreted 
as an UC Hll region, the size of the H ll region is only 2. 5 
mpc and is excited by a Bl ZAMS star (jHofner et al.lll996h . 
iFranco et all (|2000l ) quote a slightly flatter spectral energy 
distribution between 8.4 and 110 GHz, with a spectral in- 
dex of 0.95 ±0.06 which they attribute to a dens ity gradient 
prop ortional to r~^'^. Using the 1.3 cm dat a of iTesti et al.l 
l|200a ') for G9.62+0.20E. [Franco eral] (|2000l ') derives an elec- 
tron density of about 6 x 10^ cm~^ which could rise to above 
10^ cm~^ in a uniform core region. Taken together with the 
size of the Hll region, G9.62+0.20E qualifies as a hypercom- 
pact H II region, implying a very early phase of massive star 
formation. 

Component E is also well defined in a number of sul- 
phur and nitroge n bearing mole cules as well as in some or- 
ganic molecules l|Su et ahllioosi ). Near-infrared imaging of 
the star forming region shows emission at c omponents B, 
C and E l|Testi et al.lll998h . |Persi et al.1 (|2003l ) subsequently 
found that the NIR s ource near comp o nent E has the colours 
of a foreground star. lDe Buizer et al.l (J2005') find weak emis- 
sion at 11.7 /im towards the radio component E. 

In this paper we present further monitoring data on 
G9.62+0.20E for methanol masers at 6.7, 12.2, and 107 GHz. 
We present various analyses of the data and consider a pos- 
sible explanation for the variability of the masers. 



2 OBSERVATIONS AND DATA REDUCTION 

HartRAO observations: Observations at 6.7 GHz and 12.2 
GHz were made with the 26-m Hartebeesthoek telescope. 
Both receivers provide dual circular polarization. The sys- 
tem temperatures at zenith are typically 60 K at 6.7 GHz 
and 100 K at 12.2 GHz. The spectrometer provides 1024 
channels per polarization, and the bandwidths used for spec- 
troscopy of 1 MHz at 6.7 GHz and 2 MHz at 12 GHz provide 
spectral resolutions of 0.044 and 0.048 km s~^ respectively. 

For spectroscopic observing, the telescope pointing er- 
ror was determined through five short integrations at the 
cardinal half power points of the beam and on source. The 
on-source spectra were scaled up by determining the am- 
plitude correction from the pointing error using a Gaussian 
model for the main beam. 

Calibration was based on monitoring of 3C123, 3C218 
and Virgo A (w hich is bri g ht bu t partly resolved), using 
the flux scale of lOtt et al.l (|l994l ). The antenna tempera- 
ture from these continuum calibrators was measured by drift 
scans. Pointing errors in the north-south direction were mea- 
sured via drift scans at the beam half power points, and the 
on-source amplitude was corrected using the Gaussian beam 
model. The bright methanol maser source G351.42+0.64, 
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Figure 1. Time series for the continuum emission of DR-21 at 
107 GHz 

which previous monitoring has show n to exhibit only sma ll 
variations in the strongest features (jCoedhart et al.ll2004l ). 
was observed in the same way as G9. 62+0. 20 to provide 
a consistency check on the spectroscopy. These sources 
were observed frequently, usually on the same days as 
G9.62+0.20. 

ARO 12m observations: The observations with the ARO 
12m telescopsQ were done from December 3, 2003 to June 
25, 2004 to observe cycle 8 and from January 14, 2005 to 
April 22, 2005 to observe cycle 9. The peak for cycle 8 
was expected in early July 2004. Poor atmospheric condi- 
tions forced any further observations to be abandoned af- 
ter June 2004. The onset of flare 9 was expected to oc- 
cur around February 7, 2005 with a peak between March 
7-13 and to decay to its low state by middle April 2005. 
The receivers used were dual-channel SIS mixers operated in 
single-sideband mode. The backends were filter banks with 
lOOkHz and 250 KHz resolution. The observing frequency 
was 107.01385 GHz. Each observing run lasted about three 
hours, during which regular positional checking on planets 
were done. 

DR-21 was used as calibrator as well as a reference 
source and was observed in spectral line mode. Data re- 
duction was done with the CLASS package. Fig[l] shows the 
time series for DR-21. The results are for the 250 KHz res- 
olution and each point is the average over all channels. We 
estimated the flux density of DR-21 at 107 GHz u sing the 
continuum spectrum given by iRighini et all l|l976h . which 
gives a flux density of 16.7 Jy. For each observing run an 
appropriate conversion factor was calculated and applied to 
the corrected antenna temperature of G9.62+0.20E. 



3 RESULTS AND ANALYSIS 

The average or representative spectra for the three maser 
transitions are shown in Fig. [2] It is seen that, whereas 
for the 6.7 and 12.2 GHz spectra the strongest feature is 
located at about 1.3 km s~^, the corresponding feature in 

^ The Kitt Peak 12 Meter telescope is operated by the Arizona 
Radio Observatory (ARO), Steward Observatory, University of 
Arizona. 
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Figure 2. Representative spectra for the masers at 6.7, 12.2 GHz. 
The 107 GHz spectrum is an average taken over all observations. 
The features at 5.3 and 6.4 km s~^ in the 6.7 GHz spectrum are 
associated with G9.62+0.20D. 



the 107 GHz spectrum is the weaker maser feature. The 
107 GHz spectrum also shows a broad thermal component. 
G 9.62+0.20E h a s in t he pa st b een observed a t 107 GHz 
bv lVal'tts etHI l| 19951 . 1 19991 ) and lCaswell et al.l l|2000l ). The 
av erage 107 GH z spect rum in Fig. [2] resembl es the result 
of ICaswell et al.l l|2000h . The observations of IVal'tts et all 
l| 19951 . T999h show only one maser feature at -0.57 km s~^ 
and also no broad thermal component. 

The full time series for all the identifiable features in 
the 6.7 and 12.2 GHz spectra are shown in Figs. [3] and U 
For the 6.7 GHz maser, eleven different features could be 
identified in the single dish spectrum and ten for the 12.2 
GHz maser. The time series covers about 2670 days and 13 
flares. Inspection of Fig. |3] shows that while there are some 
similarities in the time series of some of the 6.7 GHz maser 
features, there also are some significant differences. 

The strongest varying 12.2 GHz masers are significantly 
more variable than the strongest varying 6.7 GHz masers. To 
quantify the variability we calculated the relative amplitude, 
defined as 
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for a number of features in the three transitions. The results 
are given in Table [1] The errors indicate the spread in the 
relative amplitudes for the last six fiares in the time series. 



No error could be given for the single 107 GHz fiare that 
has been monitored. The 12.2 GHz masers are indeed sig- 
nificantly more variable than the 6.7 and 107 GHz masers 
while the relative amplitudes of the 6.7 and 107 GHz masers 
are approximately the same. 

Figures [3] and |4] also show that the regular fiaring does 
not occur in all of the 6.7 and 12.2 GHz maser features. 
For the 12.2 GHz masers it is seen that the fiaring be- 
haviour is limited to three features namely -(-1.25 kms~^, 
+ 1.63 kms~^ and -1-2.13 kms^^. The strongest fiaring oc- 
curs at +1.25 and +1.63 km s^^ . These are also the strongest 
two features in the maser spectrum. The absence of fiaring 
behaviour outside the above velocity range is quite obvious. 
Because of the smaller variability of the 6.7 GHz masers the 
question of which features are fiaring is not as simple as for 
the 12.2 GHz masers. Simple visual inspection suggests that 
flaring occurs for the features at -0.2, +1.2, +1.8, and +2.2 
km s^^ . 

It is noteworthy that the most variable and brightest 
12.2 GHz masers, at +1.63 and +1.25 kms~^, lie close 
to each other at the northern end of a north-south chain 
of maser features as s een in the high resolution maps of 
iGoedhart et al.l l|2005bl ). As already pointed out, these are 
also the two brightest masers in the 12.2 GHz spectrum. The 
maser features that show no periodic flaring lie signiflcantly 
offset to the south from the above two features. 

We also checked for periodic variation s in the light 
curves using the Lomb-S cargle period ogram (IScarglelll98a ) 
and Davies' L-statistic (|Daviesl [l990l '). For the 12.2 GHz 
masers both methods confirm the existence of a periodic 
signal with period 244 days in the light curves of the 
+ 1.25 kms~^, +1.63kms~^ and +2.13 kms~^ features at 
a very high level of significance. For the 6.7 GHz masers 
the periodograms were significantly more complicated than 
that of the 12.2 GHz masers. A full analysis and discussion 
of the periodograms falls outside the scope of the present 
paper and will not be discussed here. Within the context of 
the present paper it is necessary to note that for the 12.2 
GHz masers, periodic (regular) fiaring is observed only in 
three features. For the 6.7 GHz masers, fiaring behaviour is 
strongest also for only a subset of the maser features, some 
of which coincide in velocity with the 12.2 GHz masers that 
show flaring behaviour. 

In Fig. [5] we show the time series for cycles 9 and 10 for 
which the 107 GHz masers have also been monitored. The 
vertical dashed lines give the position of the peak of the 
flare at 12.2 GHz for the 1.25 kms~^ feature. For cycle 9 
the scatter in the 107 GHz data is quite large and, although 
there might be a slight hint of a rise in the flux density 
towards the expected time when the flare was at its peak, no 
clear sign of a fiare can be identified. The reason for the large 
scatter in the data is most probably the poor atmospheric 
conditions that prevailed at that time and which led to the 
termination of the monitoring for cycle 9. 

Closer inspection of the time series in Figs. [3] ID and [5] 
suggests that the flare profiles for the 6.7, 12.2 and 107 GHz 
masers might be different. To investigate to what extent the 
flare profiles are similar we flrst considered the average fiare 
profiles by calculating the normalized profile 
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Table 1. Relative amplitudes (eq. [T} for selected masers at 6.7, 12.2 and 107 GHz 



6.7 GHz 12.2 GHz 107 GHz 



Velocity 
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km s~^ 
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km s~^ 


amplitude 


1.18 


0.26 ± 0.04 


1.25 


2.0 ± 0.4 


1.14 


0.31 


1.84 


0.32 ± 0.05 


1.63 


2.4 ± 0.3 






2.24 


0.33 ± 0.04 
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Figure 5. Time series for masers at 6.7, 12.2 and 107 GHz. For 
the 6.7 and 12.2 GHz masers only one velocity component has 
been considered. The solid dashed vertical lines indicate the po- 
sition of the peak of the 12.2 GHz maser. 



for the +1.25 kms^^ feature of the 12.2 GHz spectrum and 
for the corresponding feature at 6.7 GHz. For both these fea- 
tures a 30 point running average was used after detrending 
the time series and folding of the data modulo the period of 
244 days. The 107 GHz data for the one flare were scaled ac- 
cordingly. The results are shown in Fig. [6] In spite of the fact 
that the average profiles do not coincide exactly in time, it is 
remarkable to note that the profiles seem to be very similar 
for the three maser transitions. 

The same analysis as above was also carried out on in- 
dividual 6.7 and 12.2 GHz flares to examine to what extent 
individual flares have the same profile. The result is shown 
in Fig. [T] For clarity only the flares that occured after MJD 
52750 are shown and only for the features at 1.25 kms~^. 
In spite of the large errors on the 6.7 GHz data, it is seen 
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Figure 6. Comparison of the average flare profllo for the 1.24 
km feature for the 6.7, 12.2, and 107 GHz masers. For the 
6.7 and 12.2 GHz masers the flare proflle was obtained from a 30 
point running average. 
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Figure 7. Comparison of the individual normalized flare profiles 
for the 1.24 km s~^ feature for the 6.7 and 12.2 GHz masers for 
MJD > 52750 



that the flare profiles are basically the same, also for indi- 
vidual fiares. The similarity of the individual flare profiles 
suggests that the characteristics of the mechanism underly- 
ing the flaring remain the same from flare to flare and that 
it affects the 6.7 and 12.2 GHz masers in the same way. 

We finally note that within the velocity resolution of 
the observations, no velocity shifts of the different feature 
in the maser spectrum could be detected during the flares. 



4 DISCUSSION 



iGoedhart et al.l (|2005al l briefly considered a number of pos- 
sible mechanisms that might underlie the flaring behaviour 
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Figure 3. 6.7 GHz time series 



of the masers in G9.62+0.20E. These included: disturbances 
like shock waves or clumps of matter passing through the 
masing region; variations in the flux of seed photons or of 
pump photons due to stellar pulsations; periodic outbursts, 
or the effects of a binary system. While it is not possible to 
come to a final conclusion about the origin of the period- 
icity of the masers with the present data, it is possible to 



at least exclude certain mecliarusms and point to possible 
explanations for the periodicity. 

It has already been noted that the return of the masers 
to basically the same quiescent level between flares is a 
strong indication that the masing regions remain unaffected 
by whatever mechanism underlies the flaring. Mechanisms 
such as shock waves or clumps of matter passing through 
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Figure 4. 12.2 GHz time series 



the masing region are therefore excluded. This imphes that 
whatever the underlying mechanism for the periodicity, the 
coupling with the masing regions must be radiative, ie. either 
through the seed photons or through heating and cooling of 
the dust that is responsible for the pumping radiation field. 

At present the information available on the masers are 
the single dish light curves (as presented above) and the 



high re solution interferometric mapping of iGoedhart et ahl 
(|2005bl ) . We consider the light curve (flare profile) only since 
in general the light curves of periodic sources contain sig- 
nificant information about the physical processes that are 
responsible for the observed variability. We also present an 
analysis that strongly suggests that the decay part of the 
light curve might be due to the background Hll region de- 
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caying from a higher to a lower state of ionization. It is 
necessary to point out here that since the relative positions 
of the masers with respect to the H II region is not known to 
milli-arcsecond accuracy, a more complete understanding of 
all aspects of the variability of the masers in G9.62+0.20E 
cannot be given here. 

For the present discussion we will focus on the strongest 
maser features and take the average flare profile of the 12.2 
GHz masers (Fig. |6} as representative of all flares. The flare 
profile can be described by a rather steep rise followed by a 
decaying part lasting about 100 days to eventually reach a 
minimum (or possibly quiescent) state after which it appears 
to slowly increase before the next flare starts again. Since, 
as we have already argued that effects that physically affect 
the masing region can be ruled out, we are basically left 
with two possibilities to understand the light curve: stellar 
pulsation and a binary system. 

A comparison of the maser light curves as described 
above with the light curves of all classes of pulsating stars 
clearly shows that the maser light curve has a completely 
different behaviour from that of any class of pulsating stars. 
In fact, considering the physics of stellar pulsation it is hard 
to see how to produce a light curve for a pulsating star such 
that it resembles that of the masers in G9.62+0.20E. For the 
masers the light curve is strongly suggestive of a flaring state 
that is superposed on top of a base level emission. Pulsating 
stars are not in an equilibrium state and therefore do not 
have such a quiescent state on top of which pulsations or 
flares are superposed. Explaining the observed maser light 
curves as being due to pulsation of the central star would 
therefore appear to be very difficult. 

The only alternative is that of a binary system. The 
light curve does not suggest that we are dealing with an 
eclipsing effect for the masers. A class of binary systems 
which potentially can meet the requirement of providing 
radiative coupling between the source of variability and 
the masing r e gions is the colliding- wind binaries. In fact, 
IZhekov et al] {iggi) argued that, given the high binary fre- 
quency found in young stars and the observed mass loss rates 
and wind velocities, supersonically coUiding-wind systems, 
should also be found amongst pre-main sequence binary sys- 
tems. Although it is not possible with the available data to 
come to a conclusion here whether or not G9. 62-1-0. 20E in- 
deed harbours a coUiding-wind binary, it is worth exploring 



this possibility somewhat more by considering some proper- 
ties of coUiding-wind binaries. 

The main ingredient of a coUiding-wind binary is the 
two oppositely facing shocks separated by a contact dis- 
continuity. The postshock temperature is given by T^h = 
3mvl/16k (jStevens et al.lll992l : IZhekov et al.lll994l ). where 
m is the mean mass of the particles that constitute the 
wind, and u^, the wind speed. Given that the central star 
for G9.62-f 0.20E is e stimated to be o f spect ral type Bl and 
using the results of iBernabeu et al] ([198^, a wind speed 
of 800 km s~^ does not seem to be an unrealistic assump- 
tion. Using this value for v-w and assuming for simplicity a 
pure hydrog en wind, a postshock temperature of ~ 1.5 x 10^ 
K is found. iPittard et al.l l|2005h presented the emissitivity 
for hot thermal plasmas of temperatures 10^, 10 and 10* 
K, from which an extrapolation to energies below 100 eV 
suggests that the post-shock gas will emit photons from 
the visible up to X-ray energies. U sing the CHIANTI code 
(|Dere et al.ll 19971 : iLandi et aLlbOOd ). we confirmed that this 
is indeed the case. This range of photon energies obviously 
includes photons that can heat the inner edge of the cir- 
cumstellar dust as well as ionizing photons that can cause 
additional ionization in the H II region, respectively affecting 
the pumping radiation field and/or the background source 
of seed photons. Due to the lack of a numerical model it was 
not possible to calculate expected absolute fluxes of ionizing 
photons produced at the shocks as well as of photons that 
can heat the dust. A quantitative evaluation of the effects 
these photons might have on the Hll region and the dust 
was therefor not possible. 

A second property of coUiding-wind systems relevant to 
this discussion, is that the total luminosity at the shock 



scales like 



where r is the distance between the two 



stars (|Stevens et all Il992l : [Zhekov et al] Il994l '). Obviously 
this scaling property implies a modulating effect of the total 
luminosity at the shock during the orbital motion for eccen- 
tric orbits. To investigate the time dependence of , we 
used Kepler's third law to calculate the semi-major axis for 
the binary system. For this we used a mass of 17 Mq for 
the ionizing star (spectral type Bl) and arbitrarily adopted 
a mass of 8 Mq for the secondary star. Since in Kepler's 
third law the semi-major axis depends on the inverse of the 
third root of the total mass, the flnal answer is not very sen- 
sitive to the mass of the secondary. Using the period of 244 
days it then follows that the semi-major axis of the orbit is 
2.23 AU. For an elliptic orbit the radial distance between 
the two stars is given by r = a(l — e^)/(l -I- ecosd) with 
o the semi-major axis and e the eccentricity. The resulting 
orbital modulation, for eccentricities of 0.5, 0.7, and 0.9, of 
the luminosity is shown in Fig. [8] with r in AU. The position 
of periastron passage has been set at 43.5 days in all three 
cases. 

Combining the above two basic properties of coUiding- 
wind binary systems, it is clear that such systems can pro- 
vide a periodic source of photons that can heat the circum- 
stellar dust and thereby possibly affect the pumping radi- 
ation field for the masers, as well as ionizing photons that 
can cause additional ionization in the H II region. Although 
the exact projection of the masers agai nst the H II region in 
G9. 62-1-0. 20E is not known, the work o f lPhillips et"al](|l998D 
suggests that the Hll region might indeed be the source of 
seed photons for the masers. 
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Figure 9. Normalized average flare profile for the 12.2 GHz 
masers with the solid and dashed lines the fits from eq. |3] The 
vertical dashed lines indicate the time intervals which was used 
for the fitting. 



Comparison of Fig. [8]with the maser flare profile, how- 
ever, shows that, even in the case of an eccentricity of 0.5, 
the observed maser flare proflle has a decay time that is 
signiflcantly longer than what is expected simply from the 
radiation pulse produced at periastron passage. Within the 
framework of the colliding- wind binary scenario, the ob- 
served decay must then be due to either the cooling of dust 
or the recombining of the H II region (or parts thereof) from 
a higher degree of ionization to its pre-flare state. 

Usi ng the fact th at the dust cooling time is proportional 
to (|Krugelll2003l ) for the optical thin case and propor- 
tional to T~* in the optically thick case, and that in the 
optical thin case the typical cooling times for smal l grains 
with a temperature of 60 K is about 10 seconds l|Kriigell 
l2003h . an upper limit of about 10 hours for the cooling time 
is found for the optically thick case. Assuming now a grain 
temperature of lOOK and that the cooling time under opti- 
cally thin conditions is also 10 seconds, which, due to the 
dependence should actually be less, a cooling time of 
1.2 days is found for the optically thick case. Thus, even in 
the optically thick case the dust cooling time is only a small 
fraction of the decay time of the maser flare, suggesting that 
the decay of the maser flare is most probably not due to the 
cooling of dust. 

To investigate the possibility that the decaying part of 
the maser flare is due to the recombination of the H II re- 
gion or parts thereof, we note that the recombi nation time 
scale of a hydrogen plasma is given by (arie)"^ (|Osterbrockl 
I1989D . with a the recombination coefficient and Tie the elec- 
tron density. Using a = 2.95 x lO"'^^ cm"^s~^, characteristic 
decay times between 40 and 400 days are found for densi- 
ties between rie — 10® cm~^ and Ue = 10^ cm~^. This is 
signiflcantly longer than the dust cooling time and seems to 
be able to account for the observed decay time of about 100 
days. 

In view of these numbers, the question arises whether 
the decay of the maser flare can indeed be ascribed to a 
change in the seed photon flux from a recombining thermal 
plasma. In Appendix A we derive an expression for the time 
dependence of the electron density in a volume of ionized 
hydrogen going from a higher level of ionization to a lower 
equilibrium ionization state, due to the recombination of 



the plasma. It is shown that in the optically thin case, the 
intensity of the associated free-free emission decreases with 
time according to 



7^(t) oc 



1 + tanh(arze_*t) 



uq -f tanh(Q!rie_it) 



(3) 



with Uq = rie.o/jie,* > 1. n^fi is the electron density from 
where the decay starts at time t = Q, and rie,* the equilib- 
rium electron density determined by ionization balance for 
the ionizing stellar radiation. Fitting eq. |3]to the decay part 
of the flare will give an indication to what extent the decay 
of the maser flare can be interpreted as being due to the 
recombination of a hydrogen plasma. An initial non-linear 
regression analysis using the above time dependence sug- 
gested that the decay of the maser flare actually consists of 
two components that have to be fltted separately. The flrst 
component is from days 48 to 66, and the second from days 
67 to 135. 

Following the discussion in Appendix A, we thus fit- 
ted the two components with eq. |3] by systematically vary- 
ing ne,*. This allowed us to estimate «o and therefore 
ne,o ~ uoUe,*, ie. the electron density at the time when 
the ionization pulse was switched off or when its effect was 
not significant anymore. For the interval between days 48 
and 66 it was found that 7ie,o varies from 1.57 x 10^ cm~'^ 
to 2.00 X 10*^ cm"^ for rie,* ranging from 1.0 x lO'' cm"^ 
to 7.0 X 10^ cm~^. Similarly, for the interval between 
days 67 and 135, ne,o varies between 6.0 x 10^ cm~^ and 
7.3 x 10^ cm~^ for ne,* ranging from 1.0 x lO** cm"'' to 
2.0 X 10^ cm~^. For each interval the lower value for the 
range of rie,* was chosen such that ne,o is also the solu- 
tion of eq. IA5I ie. the solution for the limiting case when 
rie,* <^ riefi- The upper value for the range of n^,* was cho- 
sen to be about one third or half of the value of n^^o obtained 
for the lower value of the range of Ue,*. 

In Fig. [5] we show the fit of eq.[3]to each of the two time 
intervals of the decay of the maser flare for specific values of 
rie,^ from the abovementioned ranges of values. For days 48 
to 66 we used rie,* = 3 x 10^ cm~^ and for days 67 to 135 we 
used ne,ic = 1 x 10^ cm~^. These two values correspond to 
ne,o = 1-66 X 10® cm"'' and ne,o ~ 6.36 x 10^ cm~^, respec- 
tively. The quality of the fits is quite remarkable, especially 
given that a priori there is no information suggesting that 
the decay of the maser flare might follow the recombina- 
tion of a thermal hydrogen plasma. The values for ne,o as 
estimated from the non-linear regression are also what is 
expected for hypercompact H II regio ns. In fact, the ag ree- 
ment with the densities quoted by F ranco et al.l (|2000l ) for 
G9.62-)~0.20E is also remarkable. At the least these results 
are strongly suggestive that the decay of the maser flare 
might be due to the recombination of an H II region, or parts 
thereof, from a higher to a lower ionization state. The en- 
tire maser flare might therefore be due to changes in the 
background source. 

We also note that within the framework of the hypothe- 
sis that the maser flaring is due to changes in the background 
H II region, day 48 gives an upper limit of the time when the 
radiative pulse causing the increase in the ionization level, 
has "switched" off and the H II region is left to decay to its 
equilibrium level. Referring to Fig. |S]the start of the ffare is 
at about day 10, implying that the full width of the radia- 
tive pulse is at most about 38 days. This suggests a rather 
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sharply peaked pulse which is what is qualitatively expected 
from the colliding-wind binary scenario discussed above. 

To what extent the colliding-wind binary scenario is 
the only one that can explain the maser flaring within the 
broader framework of binary systems, is not clear yet. Our 
discussion above does not include the presence of a possi- 
ble accretion disk or even two accretion disks, as well as 
outflows. Obviously such scenarios are much more complex, 
but should nevertheless be able to explain the periodic flar- 
ing as well as the shape of the flare profiles. Given the short 
cooling time of the dust, it follows, if the entire flare pro- 
file is due to changes in the pumping radiation field, that 
the primary driving source of the variable infrared radiation 
field should basically have the same time dependence as the 
maser flare proflle. Furthermore, the physical process should 
be rather stable from orbit to orbit to explain the similarity 
of the flare profiles in the time series. Obviously, monitoring 
of G9. 62-1-0. 20E in the mid-infrared might help resolve this 
problem. 



or X-ray signal associated with G9.62+0.20E? Is the ionizing 
photon flux produced in the hot post-shock gas large enough 
that it still can have an observable effect on Hii region in 
terms of changes in the free-free emission? Are there other 
mechanisms associated with young binary systems that can 
also account for the maser flaring? It is also necessary to 
determine the exact position of the masers relative to the 
Hii region. Obviously, significantly more work still needs 
to be done before a final answer on the periodic flaring in 
G9.62-I-0.20E can be given. 
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5 SUMMARY AND CONCLUSION 

We presented the results of the monitoring of methanol 
masers in G9.62-I-0.20 at 6.7, 12.2, and 107 GHz. Like the 6.7 
and 12.2 GHz masers, the 107 GHz masers also show flar- 
ing behaviour, with a relative amplitude that is the same as 
that of the 6.7 GHz masers. It was also found that the flare 
proflle is the same for the 6.7, 12.2 and 107 GHz masers, and 
that for the 6.7 and 12.2 GHz masers the profiles of individ- 
ual flares also are the same. We have also shown that in the 
low phase of the masers, i.e. between two flaring events, the 
ratios of maser flux densities return to the same levels. From 
this behaviour we conclude that the physical conditions in 
the masing region must be relatively stable and that the 
source for the flaring behaviour lies outside the masing re- 
gion. 

Comparison of the maser light curves with that of pul- 
sating stars lead us to conclude that stellar pulsation is 
not the underlying cause for the observed maser flaring. 
This loaves a binary system as the only other option. It 
was argued that, at least qualitatively, a colliding-wind bi- 
nary can provide the mechanism for a periodic background 
source and/or a periodic pumping radiation fleld. Although 
there might be heating of the dust due to radiation from the 
shocked regions, and thus the continuum of G9.62+0.20E 
might then show variability in the IR, it seems difficult to 
explain the decay time of about 100 days as being due to 
the cooling of dust. 

It was shown that the characteristic recombination time 
of an Hii region with densities in the range 10^ cm~^ to 
10® cm~^ can explain the observed decay time of the maser 
flare. We also showed that during the intervals 48 to 66 
days and 67 to 135 days, the decay of the maser flare is 
what can be expected for the recombination of thermal plas- 
mas with densities of approximately 1.6 x 10® cm~^ and 
6.0 X 10® cm~^, respectively. These values are in very good 
agreement with densities derived independently from radio 
contirmum measurement . 

A number of questions, however, still need to be an- 
swered. Cam we find observational evidence for a colliding- 
wind binary? Is there a periodic infrared, radio continuum. 
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APPENDIX A: THE DECAY OF THE 
ELECTRON DENSITY IN AN Hii REGION 

We consider a volume element of a pure hydrogen H II region 
which is not fully ionized. The ionization rate due to the 
stellar radiation is given by F*. The total hydrogen density is 
uh while the electron density due to ionization by the stellar 
radiation field is given by ne,*, the neutral hydrogen density 
by riijo and, nn = We,* -I- n^fo. We now add an additional 
time-dependent source of ionizing radiation, described by a 



time-dependent ionization rate Tp{t). The rate equation for 
the electron density is then given by 



rfWe 2 , ,^ 



-Vp(t))nH 



(Al) 



where a is the recombination coefficient. The first term on 
the right gives the decrease of the electron density due to 
recombinations and the second term the production of elec- 
trons due to photoionizations. 

At time t = Q the additional source of ionization is 
switched off, and = ne,o > n-e,*. The equation that then 
governs the recombination of the plasma is 

dUe 



dt 



= —arig + T^nH 



(A2) 



Recombination and ionization balance for the stellar radia- 



The rate equation then becomes 

(A3) 



tion implies F^^h ~ ct 

dUe 2 , 2 

— ; — = —an. + arip ^ 
dt 

This is a Riccati equation for which analytical solutions can 
be found. For eg. IA3l the solution is given by 

uq + tanh(Qne.*t) 



1 + ito tanh(ane,*t) 



(A4) 



where uo = ne,o/ne,* > 1. It is straight forward to check 
that at t = 0, He = riefi and that for large t, ne{t) — > rie,*. 

We note that a special case of eg. IA3l is when n^,* = 0. 
The second term in eg. lA3l can then be dropped. In this case 
the solution is easily found to be 



(A5) 

4] when * ^ 



1 + ane.ot 

This solution can also be found from eq. 

Since we do not measure the electron density directly 
but rather the intensity of free-free emission from the 
plasma, we need to calculate Iv{t)- We can consider the two 
extreme cases of the Hll region being optically thick or op- 
tically thin. In the optically thick case we have 1^ ~ Bi,{T) 
which is independent of the electron density and, if T is 
constant, it follows that dl^/dt = 0. 

In the optically thin case we have 

Iu{t) = T^{t)B^{T) cx nl{t)B^{T) 

since r oc n^. 
fore that 



Iv{t) OC Tie 



For the general case (eq. 
1 -I- Uo tanh(Qne.*t) 



(A6) 

U it follows there- 



to -f laxih{ane,i,t) 
while for the special case (eq. IA5[l 

In terms of logarithms we have 

1 + Uo tanh(arze_*t) 



log/y(t) OC —2 log 
and 



Uo -l- tanh{ane,~kt) 



logI^{t) OC — 21og(l + ariefit) 



(A7) 



(A8) 



(A9) 



(AlO) 



for the general and special cases respectively. 

When doing a regression analysis using eq. IA7I it is ob- 
viously better to resort to the logarithmic representation as 
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in eg. IA9I However, when written as an equality, eg. IA9l has 
three unknown parameters: a constant, uo and rie,*- Since a 
hnear regression can only give the values of two unknowns, 
it is necessary to specify a value for one of the three un- 
knowns. In the present paper we specify the value of rie,* 
since it is known what realistic values rie,* can assume. This 
allows us then to estimate uo and therefor also rie.o- 

This paper has been typeset from a linjX/ IM^HjX file prepared 
by the author. 



